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Abstract 

The mathematical model of periodically amplitude and phase modulated si- 
nusoidal oscillation is studied, and its Fourier spectrum is given analytically. 
The Fourier spectrum of the model explains the main features of the frequency 



Oh 



I spectrum of RR Lyrae stars showing light curve modulation called the Blazhko 
O effect: among others the appearance of multiplets, the rapid decrease of their 
amplitudes in increasing orders, the asymmetry of the amplitudes of the side 
] frequency pairs, and the possibility of the occurrence of frequency doublets in- 
'— stead of triplets in the spectrum. The good agreement of the results of this 
, mathematical model with observational facts favours those physical models of 
<^ . the Blazhko effect which explain the phenomenon as a modulation of the oscil- 

■ lation with the modulation frequency, /„i. 

OO 

, -I ■ 

■ Introduction 

. In the past years contradictory descriptions have been published about the fre- 
' quency spectrum of the RR Lyrae stars with modulated light curves. Analyzing 
. the light curves of RRab stars in the large data bases of MACHO and OGLE 
^ ; (Alcock et al. 2003; Moskalik &. Poretti 2003; Collinge et al. 2006) different 
• ^ . classifications of the modulation were introduced depending on the number and 
' the separation of the modulation side-frequencies detected in the spectrum. In 
^ , variables where an equidistant triplet structure around the fundamental mode 
■ ■ frequency was found, no further multiplet structure could be identified in these 
studies. The observations focused on individual Blazhko stars, however, in- 
evitably proved the presence of higher multiplets (quintuplet, septuplet struc- 
tures) in the spectrum (Hurta et al. 2008; Jurcsik et al. 2008; Kolenberg et 
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al. 2009, Jurcsik et al. 2009a). Up to now the only mathematical model that 
aimed to describe the full spectrum with multiplets of Blazhko variables was 
published by Breger & Kolenberg (2006). However, no detailed confrontation 
with the observations of the possible predictions of this model e.g., on the 
amplitudes of the components of the multiplets has been performed. 

As the frequency spectrum implies the basic information (amplitudes and 
phase angles) on the light curve modulation, its correct knowledge and inter- 
pretation is very important for proper understanding of the Blazhko effect. If 
we find the adequate mathematical model that describes all the features of 
the frequency spectrum of the well-observed Blazhko stars (frequencies, ampli- 
tudes and phase relations), it may provide a starting-point and steady base for 
theoretical investigations. 

The frequency spectrum of periodically modulated sinusoidal oscil- 
lation 

The amplitude and phase modulated sinusoidal oscillation is given by the for- 
mula 

m{t) = a [1 -|- 6 sin {Clt + ipi)] sin [ivt + (po + c sin {fit + (^2)] , (1) 

where lo = 27r/o, /o is the fundamental frequency, O = 27r/,n, /m is the 
frequency of the modulation, a, b and c = 2Trfoq are the amplitudes of the 
oscillation, the amplitude and the phase modulations, respectively, q expresses 
the amplitude of the phase modulation relative to the fundamental period, (po, 
(fi and (p2 denote the phases of the oscillation, amplitude and phase modula- 
tions, respectively. Analytically, there is no limit for the parameters a, b and 
c{q), but because of observational constrains we confine our discussion to the 
parameter values < 5 < 1 and |c| < 7r/2. Solutions for parameters out of 
these ranges are irrelevant in asteroseismology. 

We do not consider modulation of the oscillation frequency here as the 
modulated oscillation of the form of 

m{t) = a sin[a;(l -|- rfsin(n -|- ^3))* -|- (po] = 
a sm[cot + ipo + dujt sin(Ot -|- (^3)] 

corresponds to a phase modulation with variable modulation amplitude (c = 
dhjt) and time dependent, unstable frequency spectrum. 

By suitable choice of the starting epoch, without any restriction of the 
general validity, ipo = and ^2 — ^ can be attained. If the initial epoch 
corresponds to the timing of the mid of rising branch of the phase modulation 
at the moment of the mid of the descending branch of the oscillation (taking 



B. SzeidI & J. Jurcsik 



3 



into account the reverse direction of the magnitude scale) both (pQ = and 
(p2 = fulfill. Since /o and fm can be regarded as rational numbers as the 
accuracy of their numerical value is limited by the observations, such an epoch 
should exist. We note here that with cosine representation the choice of the 
appropriate initial epoch would be more natural, it would correspond to the 
maxima of the oscillation and the phase modulation. 

Denote now ^ = ipi — ip2 the epoch independent phase difference between 
the amplitude and phase modulations. Then the time history of the modulated 
oscillation is described as follows 



m{t) = a [1 + 6 sin {^It + $)] sin {ut + c sin fit) . 



(2) 



Taking the simple trigonometric addition formula of sin(Q! + /3) and the 
Taylor-series of sin a; and cosx into account, Eq. 2 has the form 



m{t) = a (1 + 6 sin $ cos Slt + b cos $ sin Clt) 



sinwt V(-l)"^— sin2" nt + coswt V(-l)" 
' (2n)i ^ 



,2n + l 



71=0 



ra=0 



(2n + l)! 



2n+l 



nt 



(3) 



Substituting the well-known power-reduction formulae Eqs. 4 and 5 into Eq. 3 
with a = Q,t 



sin^" a 



22n-l 



e„>|:(-)'(„l".)»-- 



1 /2n 
2 



(4) 



and 



sin 



2n+l 



sin(2fc + l)a 

ft "J / 



(5) 



and then applying the simple trigonometric formulae 



sin a sin /3 = - [cos {a — (3) — cos (a -|- /?)] 
sin a cos = ^ (o; — /?) + sin (a -|- /3)] 



(6) 



cos a cos /5 = 2 (c* ~ /^) + cos (a 
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where a = ut, f3 = iflt, i = 1, 2, 3, . . . , Eq. 3 will have the form 

oo oo 

m{t) = xo sinujt + zq coswt + Xi sm{uj + iCl)t + yi sin(a; — ifl)t 



i=l i=l 



+ ^ Zi cos{uj + iVL)t + ^ Wi cos(a; — iQ)t 

i=i i=i 

(7) 

where Xi,yi, Zi and Wi {i = 1,2,3, .. .) coefficients depend only on a, b, c and 
Applying the known relations Eqs. 8 and 9 to Eq. 7: 

cos a = sin (^a + , a = {u)±iCt)t i = 0,1,2,... (8) 

and 

Xi sin ut + X2 sin (^ut + = X sin (i/t + if) , i'=LO±i^ (i = 0, 1, 2, . . . ) 

(9) 

where X'^ = X^ + and tanip = X1/X2 we arrive at the Fourier spectrum 
of the modulated sinusoidal oscillation: 



m{t) = Ao sin{Lot + Xo) + ^ A+ sin [{lo + iQ.)t + xt] 

m 

i=l 

By the outlined procedure we can derive the Aq, A'^ and A^ coefficients 
(amplitudes) as well as their phase angles, xo. xt ^""^ Xi ■ 

As an example we derive the amplitudes and phases of the triplet. In this 
case only the fc = 0, 1 terms in Eq. 4 and the fc = term in Eq. 5 should be 
considered and substituted into Eq. 3 for each n (n = 0, 1, 2, . . . ). All the other 
terms contribute to the higher members of the multiplets. In this case taking 
Eqs. 6 into account, Eq. 3 takes the form: 

m(f) = a (1 + 6 sin $ cos Slt + b cos $ sin Clt) ■ 

•■^sinwt) (-1)":-T5 — Tj^r-n - , cos2m-- 

+ coswt > (-l)"-5— — ]sinnt}+g 

> 22-(2n+l)!V n ) j ^ 

(11) 
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where g is a function of w ± kCl, k >2. 
Let the convergent series 



n=0 ^ 

and 

oo 



c" c° 

+ 



.10 



n=0 



1 + — 

4 64 2304 147456 14745600 



„io 



+ ... 

(12) 



2"n! / n+1 



8 96 3072 184320 17694720 



Since 

2n 
n 



(13) 



(2n)! 



2n 
n- 1 



(2n)! n /2n + 1 
(^n + 1' V n 



(2n+l)! 1 
(n!)2 n+1 



(n!)2 ' 
Eq. 11 takes the form: 

m(t) = a (1 + 6 sin $ cos flt + b cos $ sin ilf) • 

• [S sinujt - 2Sismujtcos2nt + {S - Si)ccoaujtsmnt] + g (14) 

If we execute the multiplications in Eq. 14 taking into account Eqs. 6, and 
disregard all the terms with u) ± kCl,k > 2, we obtain the following equation 
that describes the triplet structure: 



m{t)tnplet = a 

+ 
+ 
+ 
+ 



SsiuLot +^{S — Si)bccos^cosujt 



^56 sin $ - ^Sibsin^ ~ ^("^ ~ ^^^^ 



-Sb cos <& + is'i6 cos <& 
2 2 



^Sbsin^ - ^Sibsm^ + ^{S - Si)c 



sm{u) — 0,)t 
cos(a; - n)t (15) 
sin(a; + Q,)t 



-^Sb cos^ + ^Sibcos^ 



cos{ui + Q)t 



Application of Eqs. 8 and 9 to Eq. 15 yields the squares of the amplitudes 
as well as the phases of the triplet: 



Al = a^[S^ + ^{S - Sifb^c^ cos^ $] 



(16) 
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{S - Sifib^ + + 26c sin $) + SSib^ cos^ 4> 
{S - Sifib"^ + - 26c sin $) + 5516^ cos^ $ 
tanxo = 2 — g — 6c cos $ 



tanx+ = 



tanxi 



(5 + 5*1)6008$ 
(S'-5i)(6sin$ + c) 
(S' + 5i)6cos$ 



(17) 
(18) 
(19) 



(S'-S'i)(6sin$-c) 

Note that a remarkable relation exists for the power difference of the side 
frequencies: 

(A+f - {A^f = (5 - SifaHcsin^ (20) 

Eq. 20 shows that the asymmetry of the side frequency amplitudes depends 
basically on the phase difference between the amplitude and phase modulation 
components. 

The amplitudes and phases of the subsequent frequencies of the multiplets 
can be derived in the same way. E.g., the final result for the side frequencies 
of the quintuplet (fc = 2) is as follows. Let the convergent series 



n=0 



C2 c4 



2"n!y (n+l)(n + 2) 

c^ c8 
+ 



then 



24 384 15360 1105920 123863040 



1 



and 



si - -{S - Si - S2)bcsm^ 



+ —{S -Si- S2fb^c' + -S2{S - Si)b''^ cos^ $ 



Si + -{S - Si - S2)bcsm^ 



+ iz{S-Si- S2fb^c^ + -,S2{S - Si)b'c^ cos" $ 



1 



2„2 , 



while 



16 



tanxj 



(S-5i +5'2)bccos'^ 
45i-(5-S'i -S'2)6csin$ 
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(S-Si + S2)bccos^ 

tanx2 = 



451 + (5-51 -53)60 sin 

Note again that the power difference of the frequency pair depends simply 

on sin$: 

{A+f - {A^f = -(5 - 5i - S2)aHcsm^. 

Discussion 

In realistic cases, in Blazhko stars the harnnonics of the fundannental frequency 
are also present in the oscillation and appear in the frequency spectrum. So 
the time history of the periodically modulated oscillation is described by the 
equation 

oc 

m{t) = Ui [1 + bi sin {fit + $»)] sin [iut + Ci sin Vlt] . (21) 

i=l 

In the previous section we confined ourselves to the fundamental frequency, 
but the same deduction is valid and true for the harmonics too, simply lo should 
be replaced by ito. Therefore, we conclude that the frequency spectrum of any 
periodically modulated periodic oscillation i.e. the light curve of Blazhko stars 
is the infinite series of the multiplets of the fundamental frequency and its 
harmonics: 



00 00 



m{t) = ^ A,Q sm{iLut + x^o) + [^'^ + ^ + ^51 

+ Sill [(^'^ - J'^) ^ + ] 

We also note here that, if the sums in Eq. 22 start from z = 0, then this 
description gives a natural ground to the appearance of Aqq, and frequencies 
at j/m, which correspond to an arbitrary zero point of the scale, and the 
modulation frequency and its harmonics, respectively. Although the appearance 
of fm (and perhaps 2/„) is hardly (if at all) discernible in most of the observed 
frequency spectrum of Blazhko stars, in accurate, extended data the modulation 
frequency can be always detected with very small amplitude most probably 
because the star's physical parameters (e.g. mean luminosity, temperature) 
change very weakly during the Blazhko cycle (Jurcsik et al. 2009a, 2009b). 

A priori there is no definite connection between the parameters of the mod- 
ulation. The rough constancy of the phase differences and the systematic be- 
haviour of the amplitudes of the triplet components with increasing orders i 
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Figure 1: Variations of tlie 5 (dots), Si (circles) sums and their combination term 
(crosses) which define the amplitude decrease of the /o frequency component with 
increasing phase modulation amplitude, c. 



(see Fig. 10 in Jurcsik et al. 2009a) suggests, however, the existence of some 
connection among them. If we suppose that b, c and <& are unique parame- 
ters of the modulation in the different harmonic orders then AfJ~ and ■yt/~ 
would depend only on a^. In this case the side frequencies of the different order 
pulsation components would naturally have some common properties. 

In all probability the amplitudes of the i/o frequencies in the spectrum of 
modulated light curves run more or less similarly to that of the single periodic 
RR Lyraes. It should, however, be noted that according to Eq. 16, a^lS"^ + 
[S - SiYc^ /4\ is an upper limit (6 = 1,$ = 0) for the amplitude of /o, i.e., 
A^/a = 1 only if the modulation is pure amplitude modulation (c = 0). Any 
phase modulation component lowers the amplitude of the /o frequency. Fig. 1 
shows the dependence of 5, S\ and the upper limit of the amplitude reduction 
factor (S*^ + [S — Si)'^c? /A) on c, using /o = 2.1925 c/d oscillation frequency, 
corresponding to the pulsation frequency of DM Cyg. The phase modulation 
lowers the amplitudes of the harmonic components, too, (Aio/ai < 1 if c, > 0) 
as E.q. 16 holds for the frequencies uj = ioj, too. 

As mentioned above, the accurate and well-distributed (over both frequen- 
cies /o and fm) observations of Blazhko stars clearly show the multiplet struc- 
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ture around the fundamental mode frequency and its harmonics in their Fourier 
spectrum (Jurcsik et al. 2008, Kolenberg et al. 2009). The question arises why 
the triplet is a stril<ing feature while the higher multiplets are hardly perceptible 
in the frequency spectrum of a Blazhko star. 



c=0.1 



E 
< 




c=0.5 




c=1.0 




fo-5fm fo-'«m fo-3fm VZfm ^O'^U 'o+I'm ^O+^L 'o+Sfm 'o+^fm fo+5f„ 

frequency 



Figure 2: Amplitude decrease of the multiplet components in the different orders 

of the modulation. Artificial datasets were generated for different values of c (0.1, 
0.5, 1.0) according to Eq. 3 using a = 1, b = 0.1 and $ = 72° values and /o,/m 
frequencies of DM Cyg. The phase modulation with c = 1 amplitude corresponds 
to full amplitude of the phase modulation slightly less than 1/3 of the main period 
which is substantially larger than the amplitude of the phase modulation observed in 
Blazhko stars, which is typically about 0.1 — 0.2 pulsation phase. Even in the c = 1 
simulations the amplitudes of the higher order modulation components (/o ± kfm, 
k >3) are neglectably small. 

Although the formulae of the amplitudes of the multiplets can be exactly 
derived as it was shown in the previous section, they are fairly complicated 
expressions (besides the outlined procedure is too laborious). Their behaviour 
can be, however, easily studied through the Fourier analysis of synthetic light 



10 The frequency spectrum of periodically modulated sinusoidal oscillation 



Table 1: Amplitudes of the frequency multiplets appearing in the Fourier spectrum of 
amplitude and phase modulated sinusoidal signal according to Eq. 2 using a — l,b = 
0.1, c = 0.5 parameters for different $ values. 

^ /o ~ 5/m /o — fo — Sfm fo — ^fm fo — fm /o 

36° 0.000007 0.000134 0.002067 0.025337 0.216919 0.938674 

72° 0.000003 0.000055 0.001197 0.019320 0.196687 0.938500 

90° 0.000000 0.000032 0.001025 0.018363 0.193815 0.938470 



36° 0.000014 0.000258 0.003680 0.039023 0.273228 

72° 0.000016 0.000286 0.004054 0.042408 0.288692 
90° 0.000016 0.000289 0.004102 0.042845 0.290722 



curves. As examples, Figs. 2, 3 and Table 1 show the amplitudes of subsequent 
multiplets of the frequency spectra of Eq. 2 at different choice of the modulation 
parameters. The striking feature is that the amplitudes of the subsequent 
frequencies on both sides of the fundamental frequency approach zero rapidly. 
The degree of the decrease depends on the parameter of the phase modulation 
c, the weaker the phase modulation, the speedier the amplitude decrease is. 
Even if the phase modulation has a high value, c = 1, the amplitude of the 
fo ± 5/„i frequencies are more than three orders of magnitude less than the 
amplitude of fo (see Fig 2). For a more realistic case e.g., with c = 0.5 the 
amplitude difference between the first and fifth order modulation components 
is larger than four orders of magnitude as shown in Table 1. It is now clear that 
only the observational accuracy sets limit to the perception of higher multiplets 
in the frequency spectrum. 

One further important implication of our results is that Eq. 20 proves that 
the power difference of the side frequencies in the triplet is the physically mean- 
ingful quantity to measure the asymmetry of the triplet instead of their ampli- 
tude ratios. 

Our calculations have further serious and amazing consequences on the 
symmetry and asymmetry of the triplet components. It results also from Eq. 20 
that the triplet will only be symmetrical if any of the quantities 6, c or <f> equals 
to zero, or $ = TT. If the quantities b, cand $ have appropriate numerical values, 
one of the amplitudes of the triplet side frequencies A'^ or may become 
zero. The only solution of the {Af )"^ or (A^)"^ = second order equations 
(Eqs. 17) is 6 = c and $ = 37r/2 or $ = 7r/2. In this case the spectrum of the 
modulated oscillation is doublet instead of triplet. As an instructive example 
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Figure 3: Amplitude variations of the multiplet frequency components of peri- 
odically amplitude and phase modulated sinusoidal signal as a function of 
the phase difference between the amplitude and phase modulation components. 
^0, ^^'' j ^2 ^3 ^~ denote the amplitudes of the /o,/o ± /m,/o ± 2/^, and 
/o ± 3/m frequency components. The three panels show the results for different 
combinations of 6 and c, the amplitudes of the amplitude and phase modulations, re- 
spectively. Artificial datasets were generated using Eq. 2 with a = 1 and /o = 2.3817, 
fm = 0.06043 c/d values corresponding to the pulsation and modulation frequencies 
of DM Cyg, and Fourier analysed in order to determine the amplitude values. 



Fig. 3 shows the amplitude changes of the frequencies as a function of $ for 
a = 1, 6 = 0.1 and c = 0.1, 0.5, 1.0. In Blazhko stars, usually there is a phase 
difference between the amplitude and phase modulation components (see e.g., 
Fig 8 in Jurcsik et a! 2008, and Fig 3 in Jurcsik et al 2009a) and this explains 
why asymmetric triplets are detected in most of the cases. 

Concerning the phases, from Eq. 18 it follows that if either 6 = or c = 
the phase of the oscillation Is or tt (tanxo = 0) if the initial epoch is chosen 
as described in the beginning of Sect. 1. When the modulation is pure phase 
modulation, then tanx^ = tanx^ = 0, i.e., Xi < and xf are or tt. When 
the modulation is pure amplitude modulation, then Eq. 1 can be easily solved 
with appropriate choices of the initial epoch. If the initial epoch is chosen to 
fulfill = and ipi — 0, i.e., it corresponds to the timing of the mid of 
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rising branch of the amplitude modulation at the moment of the mid of the 
descending branch of the oscillation, the solution of Eq. 1 shows that Xi = 7r/2 
and xt — — Ti'/S- Alternatively, if the initial epoch corresponds to the maximum 
phase of the amplitude modulation i.e., ipi = n/2, then Xi — Xi — Xo fulfill. 

Eq. 20 provides another interesting result. If < <& < tt, then > A[ 
and the plot that shows the amplitude of the light variation vs. phase shift of 
maximum light during the Blazhko cycle has anticlockwise progression, while in 
the case of — tt < $ < 0, A^^ < , the progression is the opposite. In reality, 
of course, the situation is somewhat more complex because the modulations of 
the harmonic components may modify the picture. 

Conclusions 

The model of periodically modulated harmonic oscillation properly explains the 
main features of the frequency spectrum of Blazhko RR Lyrae stars. It pre- 
dicts the infinite multiplet systems around the fundamental frequency and its 
harmonics, accounts for the often detected asymmetry of the amplitudes of the 
side frequency pairs in the triplets (quintuplets) and for the rapid decrease of 
the amplitudes of the multiplets with increasing orders. 

Our results also indicate that the suggested classification schemes of Blazhko 
stars based on their frequency spectrum (Alcock et al. 2003, Moskalik and 
Poretti 2003) are dubious. For example, the frequency doublets, which are 
generally interpreted with a nonradial frequency component close to the radial 
frequency, can naturally take its origin from amplitude and phase modulations 
of a single oscillation. Our analysis also show that the asymmetry of the triplets 
are a natural consequence of the mixture of amplitude and phase modulations. 
We thus conclude that the occurrence of multiplets does not necessarily imply 
the presence of more than a single oscillation, and prefers those physical models 
which connect the phenomenon to one modulation frequency, /„. Recently, 
Stothers (2006) proposed such a model for Blazhko RR Lyrae stars. In his 
model the cyclic changes in the strength of the envelope convection gives rise 
to the modulation of the periodic oscillation of these stars. 

In the near future we plan further investigations to exploit the potentials of 
this model and to find connections between the amplitudes and phase angles in 
the frequency spectrum. We also plan to carry out a more profound comparison 
with observational facts. 
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' quency spectrum of the RR Lyrae stars with modulated light curves. Analyzing 
. the light curves of RRab stars in the large data bases of MACHO and OGLE 
^ ; (Alcock et al. 2003; Moskalik &. Poretti 2003; Collinge et al. 2006) different 
• ^ . classifications of the modulation were introduced depending on the number and 
' the separation of the modulation side-frequencies detected in the spectrum. In 
^ , variables where an equidistant triplet structure around the fundamental mode 
■ ■ frequency was found, no further multiplet structure could be identified in these 
studies. The observations focused on individual Blazhko stars, however, in- 
evitably proved the presence of higher multiplets (quintuplet, septuplet struc- 
tures) in the spectrum (Hurta et al. 2008; Jurcsik et al. 2008; Kolenberg et 
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al. 2009, Jurcsik et al. 2009a). Up to now the only mathematical model that 
aimed to describe the full spectrum with multiplets of Blazhko variables was 
published by Breger & Kolenberg (2006). However, no detailed confrontation 
with the observations of the possible predictions of this model e.g., on the 
amplitudes of the components of the multiplets has been performed. 

As the frequency spectrum implies the basic information (amplitudes and 
phase angles) on the light curve modulation, its correct knowledge and inter- 
pretation is very important for proper understanding of the Blazhko effect. If 
we find the adequate mathematical model that describes all the features of 
the frequency spectrum of the well-observed Blazhko stars (frequencies, ampli- 
tudes and phase relations), it may provide a starting-point and steady base for 
theoretical investigations. 

The frequency spectrum of periodically modulated sinusoidal oscil- 
lation 

The amplitude and phase modulated sinusoidal oscillation is given by the for- 
mula 

m{t) = a [1 -|- 6 sin {Clt + ipi)] sin [ivt + (po + c sin {fit + (^2)] , (1) 

where lo = 27r/o, /o is the fundamental frequency, O = 27r/,n, /m is the 
frequency of the modulation, a, b and c = 2Trfoq are the amplitudes of the 
oscillation, the amplitude and the phase modulations, respectively, q expresses 
the amplitude of the phase modulation relative to the fundamental period, (po, 
(fi and (p2 denote the phases of the oscillation, amplitude and phase modula- 
tions, respectively. Analytically, there is no limit for the parameters a, b and 
c{q), but because of observational constrains we confine our discussion to the 
parameter values < 5 < 1 and |c| < 7r/2. Solutions for parameters out of 
these ranges are irrelevant in asteroseismology. 

We do not consider modulation of the oscillation frequency here as the 
modulated oscillation of the form of 

m{t) = a sin[a;(l -|- rfsin(n -|- ^3))* -|- (po] = 
a sm[cot + ipo + dujt sin(Ot -|- (^3)] 

corresponds to a phase modulation with variable modulation amplitude (c = 
dhjt) and time dependent, unstable frequency spectrum. 

By suitable choice of the starting epoch, without any restriction of the 
general validity, ipo = and ^2 — ^ can be attained. If the initial epoch 
corresponds to the timing of the mid of rising branch of the phase modulation 
at the moment of the mid of the descending branch of the oscillation (taking 
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into account the reverse direction of the magnitude scale) both (pQ = and 
(p2 = fulfill. Since /o and fm can be regarded as rational numbers as the 
accuracy of their numerical value is limited by the observations, such an epoch 
should exist. We note here that with cosine representation the choice of the 
appropriate initial epoch would be more natural, it would correspond to the 
maxima of the oscillation and the phase modulation. 

Denote now ^ = ipi — ip2 the epoch independent phase difference between 
the amplitude and phase modulations. Then the time history of the modulated 
oscillation is described as follows 



m{t) = a [1 + 6 sin {^It + $)] sin {ut + c sin fit) . 



(2) 



Taking the simple trigonometric addition formula of sin(Q! + /3) and the 
Taylor-series of sin a; and cosx into account, Eq. 2 has the form 



m{t) = a (1 + 6 sin $ cos Slt + b cos $ sin Clt) 



sinwt V(-l)"^— sin2" nt + coswt V(-l)" 
' (2n)i ^ 



,2n + l 



71=0 



ra=0 



(2n + l)! 



2n+l 



nt 



(3) 



Substituting the well-known power-reduction formulae Eqs. 4 and 5 into Eq. 3 
with a = Q,t 



sin^" a 



22n-l 



e„>|:(-)'(„l".)»-- 



1 /2n 
2 



(4) 



and 



sin 



2n+l 



sin(2fc + l)a 

ft "J / 



(5) 



and then applying the simple trigonometric formulae 



sin a sin /3 = - [cos {a — (3) — cos (a -|- /?)] 
sin a cos = ^ (o; — /?) + sin (a -|- /3)] 



(6) 



cos a cos /5 = 2 (c* ~ /^) + cos (a 



4 
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where a = ut, f3 = iflt, i = 1, 2, 3, . . . , Eq. 3 will have the form 

oo oo 

m{t) = xo sinujt + zq coswt + Xi sm{uj + iCl)t + yi sin(a; — ifl)t 



i=l i=l 



+ ^ Zi cos{uj + iVL)t + ^ Wi cos(a; — iQ)t 

i=i i=i 

(7) 

where Xi,yi, Zi and Wi {i = 1,2,3, .. .) coefficients depend only on a, b, c and 
Applying the known relations Eqs. 8 and 9 to Eq. 7: 

cos a = sin (^a + , a = {u)±iCt)t i = 0,1,2,... (8) 

and 

Xi sin ut + X2 sin (^ut + = X sin (i/t + if) , i'=LO±i^ (i = 0, 1, 2, . . . ) 

(9) 

where X'^ = X^ + and tanip = X1/X2 we arrive at the Fourier spectrum 
of the modulated sinusoidal oscillation: 



m{t) = Ao sin{Lot + Xo) + ^ A+ sin [{lo + iQ.)t + xt] 

m 

i=l 

By the outlined procedure we can derive the Aq, A'^ and A^ coefficients 
(amplitudes) as well as their phase angles, xo. xt ^""^ Xi ■ 

As an example we derive the amplitudes and phases of the triplet. In this 
case only the fc = 0, 1 terms in Eq. 4 and the fc = term in Eq. 5 should be 
considered and substituted into Eq. 3 for each n (n = 0, 1, 2, . . . ). All the other 
terms contribute to the higher members of the multiplets. In this case taking 
Eqs. 6 into account, Eq. 3 takes the form: 

m(f) = a (1 + 6 sin $ cos Slt + b cos $ sin Clt) ■ 

•■^sinwt) (-1)":-T5 — Tj^r-n - , cos2m-- 

+ coswt > (-l)"-5— — ]sinnt}+g 

> 22-(2n+l)!V n ) j ^ 

(11) 
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where g is a function of w ± kCl, k >2. 
Let the convergent series 



n=0 ^ 

and 

oo 



c" c° 

+ 



.10 



n=0 



1 + — 

4 64 2304 147456 14745600 



„io 



+ ... 

(12) 



2"n! / n+1 



8 96 3072 184320 17694720 



Since 

2n 
n 



(13) 



(2n)! 



2n 
n- 1 



(2n)! n /2n + 1 
(^n + 1' V n 



(2n+l)! 1 
(n!)2 n+1 



(n!)2 ' 
Eq. 11 takes the form: 

m(t) = a (1 + 6 sin $ cos flt + b cos $ sin ilf) • 

• [S sinujt - 2Sismujtcos2nt + {S - Si)ccoaujtsmnt] + g (14) 

If we execute the multiplications in Eq. 14 taking into account Eqs. 6, and 
disregard all the terms with u) ± kCl,k > 2, we obtain the following equation 
that describes the triplet structure: 



m{t)tnplet = a 

+ 
+ 
+ 
+ 



SsiuLot +^{S — Si)bccos^cosujt 



^56 sin $ - ^Sibsin^ ~ ^("^ ~ ^^^^ 



-Sb cos <& + is'i6 cos <& 
2 2 



^Sbsin^ - ^Sibsm^ + ^{S - Si)c 



sm{u) — 0,)t 
cos(a; - n)t (15) 
sin(a; + Q,)t 



-^Sb cos^ + ^Sibcos^ 



cos{ui + Q)t 



Application of Eqs. 8 and 9 to Eq. 15 yields the squares of the amplitudes 
as well as the phases of the triplet: 



Al = a^[S^ + ^{S - Sifb^c^ cos^ $] 



(16) 



6 
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{S - Sifib^ + + 26c sin $) + SSib^ cos^ 4> 
{S - Sifib"^ + - 26c sin $) + 5516^ cos^ $ 
tanxo = 2 — g — 6c cos $ 



tanx+ = 



tanxi 



(5 + 5*1)6008$ 
(S'-5i)(6sin$ + c) 
(S' + 5i)6cos$ 



(17) 
(18) 
(19) 



(S'-S'i)(6sin$-c) 

Note that a remarkable relation exists for the power difference of the side 
frequencies: 

(A+f - {A^f = (5 - SifaHcsin^ (20) 

Eq. 20 shows that the asymmetry of the side frequency amplitudes depends 
basically on the phase difference between the amplitude and phase modulation 
components. 

The amplitudes and phases of the subsequent frequencies of the multiplets 
can be derived in the same way. E.g., the final result for the side frequencies 
of the quintuplet (fc = 2) is as follows. Let the convergent series 



n=0 



C2 c4 



2"n!y (n+l)(n + 2) 

c^ c8 
+ 



then 



24 384 15360 1105920 123863040 



1 



and 



si - -{S - Si - S2)bcsm^ 



+ —{S -Si- S2fb^c' + -S2{S - Si)b''^ cos^ $ 



Si + -{S - Si - S2)bcsm^ 



+ iz{S-Si- S2fb^c^ + -,S2{S - Si)b'c^ cos" $ 



1 



2„2 , 



while 



16 



tanxj 



(S-5i +5'2)bccos'^ 
45i-(5-S'i -S'2)6csin$ 
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(S-Si + S2)bccos^ 

tanx2 = 



451 + (5-51 -53)60 sin 

Note again that the power difference of the frequency pair depends simply 

on sin$: 

{A+f - {A^f = -(5 - 5i - S2)aHcsm^. 

Discussion 

In realistic cases, in Blazhko stars the harnnonics of the fundannental frequency 
are also present in the oscillation and appear in the frequency spectrum. So 
the time history of the periodically modulated oscillation is described by the 
equation 

oc 

m{t) = Ui [1 + bi sin {fit + $»)] sin [iut + Ci sin Vlt] . (21) 

i=l 

In the previous section we confined ourselves to the fundamental frequency, 
but the same deduction is valid and true for the harmonics too, simply lo should 
be replaced by ito. Therefore, we conclude that the frequency spectrum of any 
periodically modulated periodic oscillation i.e. the light curve of Blazhko stars 
is the infinite series of the multiplets of the fundamental frequency and its 
harmonics: 



00 00 



m{t) = ^ A,Q sm{iLut + x^o) + [^'^ + ^ + ^51 

+ Sill [(^'^ - J'^) ^ + ] 

We also note here that, if the sums in Eq. 22 start from z = 0, then this 
description gives a natural ground to the appearance of Aqq, and frequencies 
at j/m, which correspond to an arbitrary zero point of the scale, and the 
modulation frequency and its harmonics, respectively. Although the appearance 
of fm (and perhaps 2/„) is hardly (if at all) discernible in most of the observed 
frequency spectrum of Blazhko stars, in accurate, extended data the modulation 
frequency can be always detected with very small amplitude most probably 
because the star's physical parameters (e.g. mean luminosity, temperature) 
change very weakly during the Blazhko cycle (Jurcsik et al. 2009a, 2009b). 

A priori there is no definite connection between the parameters of the mod- 
ulation. The rough constancy of the phase differences and the systematic be- 
haviour of the amplitudes of the triplet components with increasing orders i 
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Figure 1: Variations of tlie 5 (dots), Si (circles) sums and their combination term 
(crosses) which define the amplitude decrease of the /o frequency component with 
increasing phase modulation amplitude, c. 



(see Fig. 10 in Jurcsik et al. 2009a) suggests, however, the existence of some 
connection among them. If we suppose that b, c and <& are unique parame- 
ters of the modulation in the different harmonic orders then AfJ~ and ■yt/~ 
would depend only on a^. In this case the side frequencies of the different order 
pulsation components would naturally have some common properties. 

In all probability the amplitudes of the i/o frequencies in the spectrum of 
modulated light curves run more or less similarly to that of the single periodic 
RR Lyraes. It should, however, be noted that according to Eq. 16, a^lS"^ + 
[S - SiYc^ /4\ is an upper limit (6 = 1,$ = 0) for the amplitude of /o, i.e., 
A^/a = 1 only if the modulation is pure amplitude modulation (c = 0). Any 
phase modulation component lowers the amplitude of the /o frequency. Fig. 1 
shows the dependence of 5, S\ and the upper limit of the amplitude reduction 
factor (S*^ + [S — Si)'^c? /A) on c, using /o = 2.1925 c/d oscillation frequency, 
corresponding to the pulsation frequency of DM Cyg. The phase modulation 
lowers the amplitudes of the harmonic components, too, (Aio/ai < 1 if c, > 0) 
as E.q. 16 holds for the frequencies uj = ioj, too. 

As mentioned above, the accurate and well-distributed (over both frequen- 
cies /o and fm) observations of Blazhko stars clearly show the multiplet struc- 
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ture around the fundamental mode frequency and its harmonics in their Fourier 
spectrum (Jurcsik et al. 2008, Kolenberg et al. 2009). The question arises why 
the triplet is a stril<ing feature while the higher multiplets are hardly perceptible 
in the frequency spectrum of a Blazhko star. 



c=0.1 



E 
< 




c=0.5 




c=1.0 




fo-5fm fo-'«m fo-3fm VZfm ^O'^U 'o+I'm ^O+^L 'o+Sfm 'o+^fm fo+5f„ 

frequency 



Figure 2: Amplitude decrease of the multiplet components in the different orders 

of the modulation. Artificial datasets were generated for different values of c (0.1, 
0.5, 1.0) according to Eq. 3 using a = 1, b = 0.1 and $ = 72° values and /o,/m 
frequencies of DM Cyg. The phase modulation with c = 1 amplitude corresponds 
to full amplitude of the phase modulation slightly less than 1/3 of the main period 
which is substantially larger than the amplitude of the phase modulation observed in 
Blazhko stars, which is typically about 0.1 — 0.2 pulsation phase. Even in the c = 1 
simulations the amplitudes of the higher order modulation components (/o ± kfm, 
k >3) are neglectably small. 

Although the formulae of the amplitudes of the multiplets can be exactly 
derived as it was shown in the previous section, they are fairly complicated 
expressions (besides the outlined procedure is too laborious). Their behaviour 
can be, however, easily studied through the Fourier analysis of synthetic light 
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Table 1: Amplitudes of the frequency multiplets appearing in the Fourier spectrum of 
amplitude and phase modulated sinusoidal signal according to Eq. 2 using a — l,b = 
0.1, c = 0.5 parameters for different $ values. 

^ /o ~ 5/m /o — fo — Sfm fo — ^fm fo — fm /o 

36° 0.000007 0.000134 0.002067 0.025337 0.216919 0.938674 

72° 0.000003 0.000055 0.001197 0.019320 0.196687 0.938500 

90° 0.000000 0.000032 0.001025 0.018363 0.193815 0.938470 



36° 0.000014 0.000258 0.003680 0.039023 0.273228 

72° 0.000016 0.000286 0.004054 0.042408 0.288692 
90° 0.000016 0.000289 0.004102 0.042845 0.290722 



curves. As examples, Figs. 2, 3 and Table 1 show the amplitudes of subsequent 
multiplets of the frequency spectra of Eq. 2 at different choice of the modulation 
parameters. The striking feature is that the amplitudes of the subsequent 
frequencies on both sides of the fundamental frequency approach zero rapidly. 
The degree of the decrease depends on the parameter of the phase modulation 
c, the weaker the phase modulation, the speedier the amplitude decrease is. 
Even if the phase modulation has a high value, c = 1, the amplitude of the 
fo ± 5/„i frequencies are more than three orders of magnitude less than the 
amplitude of fo (see Fig 2). For a more realistic case e.g., with c = 0.5 the 
amplitude difference between the first and fifth order modulation components 
is larger than four orders of magnitude as shown in Table 1. It is now clear that 
only the observational accuracy sets limit to the perception of higher multiplets 
in the frequency spectrum. 

One further important implication of our results is that Eq. 20 proves that 
the power difference of the side frequencies in the triplet is the physically mean- 
ingful quantity to measure the asymmetry of the triplet instead of their ampli- 
tude ratios. 

Our calculations have further serious and amazing consequences on the 
symmetry and asymmetry of the triplet components. It results also from Eq. 20 
that the triplet will only be symmetrical if any of the quantities 6, c or <f> equals 
to zero, or $ = TT. If the quantities b, cand $ have appropriate numerical values, 
one of the amplitudes of the triplet side frequencies A'^ or may become 
zero. The only solution of the {Af )"^ or (A^)"^ = second order equations 
(Eqs. 17) is 6 = c and $ = 37r/2 or $ = 7r/2. In this case the spectrum of the 
modulated oscillation is doublet instead of triplet. As an instructive example 
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Figure 3: Amplitude variations of the multiplet frequency components of peri- 
odically amplitude and phase modulated sinusoidal signal as a function of 
the phase difference between the amplitude and phase modulation components. 
^0, ^^'' j ^2 ^3 ^~ denote the amplitudes of the /o,/o ± /m,/o ± 2/^, and 
/o ± 3/m frequency components. The three panels show the results for different 
combinations of 6 and c, the amplitudes of the amplitude and phase modulations, re- 
spectively. Artificial datasets were generated using Eq. 2 with a = 1 and /o = 2.3817, 
fm = 0.06043 c/d values corresponding to the pulsation and modulation frequencies 
of DM Cyg, and Fourier analysed in order to determine the amplitude values. 



Fig. 3 shows the amplitude changes of the frequencies as a function of $ for 
a = 1, 6 = 0.1 and c = 0.1, 0.5, 1.0. In Blazhko stars, usually there is a phase 
difference between the amplitude and phase modulation components (see e.g., 
Fig 8 in Jurcsik et a! 2008, and Fig 3 in Jurcsik et al 2009a) and this explains 
why asymmetric triplets are detected in most of the cases. 

Concerning the phases, from Eq. 18 it follows that if either 6 = or c = 
the phase of the oscillation Is or tt (tanxo = 0) if the initial epoch is chosen 
as described in the beginning of Sect. 1. When the modulation is pure phase 
modulation, then tanx^ = tanx^ = 0, i.e., Xi < and xf are or tt. When 
the modulation is pure amplitude modulation, then Eq. 1 can be easily solved 
with appropriate choices of the initial epoch. If the initial epoch is chosen to 
fulfill = and ipi — 0, i.e., it corresponds to the timing of the mid of 
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rising branch of the amplitude modulation at the moment of the mid of the 
descending branch of the oscillation, the solution of Eq. 1 shows that Xi = 7r/2 
and xt — — Ti'/S- Alternatively, if the initial epoch corresponds to the maximum 
phase of the amplitude modulation i.e., ipi = n/2, then Xi — Xi — Xo fulfill. 

Eq. 20 provides another interesting result. If < <& < tt, then > A[ 
and the plot that shows the amplitude of the light variation vs. phase shift of 
maximum light during the Blazhko cycle has anticlockwise progression, while in 
the case of — tt < $ < 0, A^^ < , the progression is the opposite. In reality, 
of course, the situation is somewhat more complex because the modulations of 
the harmonic components may modify the picture. 

Conclusions 

The model of periodically modulated harmonic oscillation properly explains the 
main features of the frequency spectrum of Blazhko RR Lyrae stars. It pre- 
dicts the infinite multiplet systems around the fundamental frequency and its 
harmonics, accounts for the often detected asymmetry of the amplitudes of the 
side frequency pairs in the triplets (quintuplets) and for the rapid decrease of 
the amplitudes of the multiplets with increasing orders. 

Our results also indicate that the suggested classification schemes of Blazhko 
stars based on their frequency spectrum (Alcock et al. 2003, Moskalik and 
Poretti 2003) are dubious. For example, the frequency doublets, which are 
generally interpreted with a nonradial frequency component close to the radial 
frequency, can naturally take its origin from amplitude and phase modulations 
of a single oscillation. Our analysis also show that the asymmetry of the triplets 
are a natural consequence of the mixture of amplitude and phase modulations. 
We thus conclude that the occurrence of multiplets does not necessarily imply 
the presence of more than a single oscillation, and prefers those physical models 
which connect the phenomenon to one modulation frequency, /„. Recently, 
Stothers (2006) proposed such a model for Blazhko RR Lyrae stars. In his 
model the cyclic changes in the strength of the envelope convection gives rise 
to the modulation of the periodic oscillation of these stars. 

In the near future we plan further investigations to exploit the potentials of 
this model and to find connections between the amplitudes and phase angles in 
the frequency spectrum. We also plan to carry out a more profound comparison 
with observational facts. 
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